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The adsorption characteristics of alkali, alkaline earth and transition metal adatoms on silicene, a 
graphene-like monolayer structure of silicon, are analyzed by means of first-principles calculations. 
In contrast to graphene, interaction between the metal atoms and the silicene surface is quite strong 
due to its highly reactive buckled hexagonal structure. In addition to structural properties, we 
also calculate the electronic band dispersion, net magnetic moment, charge transfer, workfunction 
and dipole moment of the metal adsorbed silicene sheets. Alkali metals, Li, Na and K, adsorb 
to hollow site without any lattice distortion. As a consequence of the significant charge transfer 
from alkalis to silicene metalization of silicene takes place. Trends directly related to atomic size, 
adsorption height, workfunction and dipole moment of the silicene/alkali adatom system are also 
revealed. We found that the adsorption of alkaline earth metals on silicene are entirely different 
from their adsorption on graphene. The adsorption of Be, Mg and Ca turns silicene into a narrow 
gap semiconductor. Adsorption characteristics of eight transition metals Ti, V, Cr, Mn, Fe, Co, Mo 
and W are also investigated. As a result of their partially occupied d orbital, transition metals show 
diverse structural, electronic and magnetic properties. Upon the adsorption of transition metals, 
depending on the adatom type and atomic radius, the system can exhibit metal, half-metal and 
semiconducting behavior. For all metal adsorbates the direction of the charge transfer is from 
adsorbate to silicene, because of its high surface reactivity. Our results indicate that the reactive 
crystal structure of silicene provides a rich playground for functionalization at nanoscale. 

PACS numbers: 81.16.Pr, 68.65.Pq, 66. 30. Pa, 81.05.ue 



I. INTRODUCTION 

Recent advances in controllable synthesis and charac- 
terization of nanoscale materials, have opened up im- 
portant possibilities for the investigation of ultra-thin 
two-dimensional systems. Chiefly the research efforts di- 
rected towards graphene^ have dominated the new era 
of two-dimensional materials. Many exceptional features 
of atomically thin graphene layers such as massless Dirac 
fermions, strength of the lattice structure, high ther- 
mal conductivity and half-integer Hall conductance have 
been revealed so farPHU J n spite of its unique proper- 
ties, due to the lack of a band gap and its weak light 
adsorption, graphene research efforts have focused on 
graphene composites over the past five years. Studies 
have demonstrated the existence of several chemically 
converted graphene structures such as grapheneoxide 
(GO)EH2lgraphane (CH)CSHIllfluorographene (CF)J±5HH] 
and chlorographene (CCl)P^^ High- quality insulating 
behavior, thermal stability and extraordinary mechani- 
cal strength of fluorographene (CF) has inspired intense 
research on halogenated graphene derivatives. 

Unusual properties of graphene promising for a variety 
of novel applicationsS2H2S have also triggered significant 
interest in one or several atom-thick honeycomb struc- 
tures of binary compounds. Early experimental studies 
aiming to synthesize and characterize novel monolayer 
materials have revealed that graphene-like sheets of BN 
are also stable pSES Though BN has the same planar 
structure as graphene due to the ionic character of B- 
N bonds, BN crystal is a wide band gap insulator with 
an energy gap of 4.6 eVP^HUl The perfect lattice match- 



ing between graphene and BN layers make it possible 
to construct nanoscale devices.^ Following the synthesis 
of hexagonal monolayer of ZnOp^the II- VI metal-oxide 
analogue of graphene, it was also predicted that ZnO 
nanoribbons have ferromagnetic order in their ground 
statefSl In addition to these, it was reported that com- 
mon solvents can be used to exfoliate transition metal 
dichalcogenides and oxides such as M0S2, WS2, MoSe2, 
MoTe 2 , TaSe 2) NbSe^ NiTe 2l Bi 2 Te 3 and NbSe 2 to ob- 
tain single layers ESHIIl Most recently, possibility of var- 
ious combinations of MX 2 type single-layer transition- 
metal oxides and dichalcogenides, stable even in free- 
standing form, was also predicted!^! 

The recent synthesis of silicene JSH2H the silicon 
analogue of graphene has opened a new avenue to 
nanoscale material research. Though the nanotube^ 
and fullerene 48 forms of silicon were synthesized earlier, 
monolayer silicon was presumed not to exist in a free- 
standing form. Early theoretical works pointed out that 
silicone is a scmimctal with linearly crossing bands and 
it has a buckled crystal structure that stems from sp 3 
hybridizationP2Ml Similar to graphene, the hexagonal 
lattice symmetry of silicene exhibits a pair of inequiv- 
alent valleys in the vicinity of the vertices of K and 
K' symmetry points. Moreover, the experimental real- 
ization of the transformation of thin films of wurtzite 
(WZ) materials into the graphene-like thin film struc- 
ture is another evidence for the existence of monolayer 
structures of Si and Ge.^3 Recent theoretical studies have 
revealed several remarkable features of silicene such as a 
large spin-orbit gap at the Dirac point, 52 experimentally 
accessible quantum spin Hall effect,^ transition from a 
topological insulating phase to a band insulator that can 



2 




VI K r M 



FIG. 1: (Color online) (a) Top view of the honeycomb lattice 
of silicene and it's unitcell shown by black parallelogram, (b) 
Tilted view of 2x2 supercell of silicene. Structural parameters: 
buckling (5), Si-Si-Si angle and Si-Si distance are indicated, 
(c) Electronic band dispersion of silicene (with and without 
SOC) and band decomposed charge densities of VB and CB 
at r, M and K symmetry points. 

be induced by an electric fielcP^ and electrically tunable 
band gapP^ In addition to unique insulator phases such 
as quantum spin Hall, quantum anomalous Hall and band 
insulator phases, the emergence of valley-polarized metal 
phase was also reported very recently. 56 It appears that 
silicene will be a possible graphene replacement not only 
due to its graphene-like features but also because of its 
compatibility to existing silicon-based electronic devices. 

In this paper, motivated by the very recent experi- 
mental realizations of monolayer silicene we inves- 
tigate how alkali, alkaline earth, and transition metal 
atoms interact with monolayer freestanding silicene. De- 
tails of computational methodology are described in Sec. 
II. Characteristic properties of monolayer silicene and 
graphene are compared briefly. Our results on structural 
and electronic properties of metal adatom adsorbed sil- 
icene are presented in Sec. IV. Conclusions and a sum- 
mary of our results are given in Sec. V. 



II. COMPUTATIONAL METHODOLOGY 

To investigate the adsorption characteristics of al- 
kali metals and transition metals on a monolayer hon- 



TABLE I: Calculated values for graphene, and silicene. These 
are lattice constant (a); Si-Si (or C-C) bond distance (d); 
thickness of the layer (t); workfunction ($); cohesive energy 
per unit cell (E co h); in-plane stiffness (C); optical phonon 
modes at the F point. 
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eycomb structure of silicone we employ first-principles 
calculations^ using the projector augmented wave 
(PAW) method 5 ® implemented in VASP code. Electronic 
exchange-correlation effects are simulated using the spin 
polarized local density approximation^ (LDA). For the 
plane-wave basis set, the kinetic energy cutoff is taken to 
be h 2 \k + G\ 2 /2m = 500 eV. Brillouin zone (BZ) sam- 
pling is determined after extensive convergence analysis. 
In the self-consistent potential, total energy and binding 
energy calculations with a (6x6x1) supercell of silicene 
sheet a set of (5x5x1) k-point sampling is used for BZ 
integration. For partial occupancies the Gaussian smear- 
ing method is used. The convergence criterion of our 
self consistent calculations for ionic relaxations is 10~ 5 
eV between two consecutive steps. By using the conju- 
gate gradient method, all atomic positions and the size of 
the unitcell were optimized until the atomic forces were 
less than 0.05 eV/A. Pressures on the lattice unit cell are 
decreased to values less than 1 kB. Adatom adsorbed sil- 
icene monolayers are treated using a supercell geometry, 
where a minimum of 15 A vacuum spacing is kept be- 
tween the adjacent silicene layers. Diffusion pathway of 
adatoms are calculated for 10 different adsorption points 
on (4x4x1) silicene supercell. 

The cohesive energy of silicene (also for graphene) per 
unit cell relative to free Si atom, given in Table IT] is ob- 
tained from E coh = 2Ef - E% iHoene , where Eip is the 
total energy of single free Si and E^ tllcene is the total en- 
ergy of silicene. Here total energy of single atoms are cal- 
culated by considering their magnetic ground state. As 
for the adsorption energy of a metal adatom, one can use 
the formula E Ads = E Ad + E^ licene - E^ iUcene+Ad where 

gSilicene^ E Ad and fiSMcene+Ad &Te thg total energies Q f 

(6x6x1) supercell of silicene, isolated single adatom and 
silicene+adatom system, respectively. 

Most of the adatom adsorption result in a net 
electrical-dipolc moment perpendicular to the plane. 
Therefore ground state electronic structure, magnetic 
state and workfunction are calculated by applying a 
dipole correctional to eliminate the artificial electrostatic 
field between the periodic supercells. To obtain continu- 
ous density of states curves and to determine the energy 
bandgap (E g ) smearing with 0.2 and 0.001 eV is used, re- 
spectively. For the charge transfer analysis, the effective 
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FIG. 2: (Color online) Preferable adsorption sites: Hollow, 
Top, Hill and Bridge on a silicene lattice. 

charge on the atoms are obtained by the Bader method.^ 

III. SILICENE 

Though crystalline silicon has the diamond structure 
and no layered form exists in nature, very recent experi- 
mental studies have reported the successful synthesis of a 
monolayer of silicon, called silicene, by the application of 
various deposition techniques. Similar to graphene, sil- 
icene can be viewed as a bipartite lattice composed of two 
inter-penetrating triangular sublattices of silicon atoms. 
Since it bonds between silicon atoms are weaker than in 
the case of the carbon atoms, planarity is destabilized and 
therefore silicone atoms are buckled in a silicene crystal. 
As shown in Fig. [ljb) the buckling (perpendicular dis- 
tance between these two Si planes) is 0.49 A. Upon the 
formation of sp 3 bonded honeycomb lattice, the covalent 
bond length of Si-Si is 2.25 A. 

Two dimensional silicone sheet is a semimetal because 
the valence and conduction bands touch at the Fermi 
level. It was predicted earlier that similar to graphene, 
silicene has also linearly crossing bands at the K (and K 1 ) 
symmetry points and charge carriers in graphene behave 
like relativistic particles with a conical energy spec trum 
with Fermi velocity Vp = 10 6 m/s like in graphene ! 49 l 5 ° l 
In Fig. [ljc), the electronic band structure of perfect sil- 
icene is presented. Linear n and n* bands that cross 
at the K symmetry point are responsible for the exis- 
tence of massless Dirac fermions in silicene. Due to the 
degeneracy of the valence band (VB) maxima and the 
conduction band (CB) minima at the K point, corre- 
sponding states have the same ionization potential and 
electron affinity. Therefore one can expect the observa- 
tion of similar unique properties of graphene in silicene. 
The calculated energy band gaps at M and V symmetry 
points are 1.64 and 3.29 eV, respectively. At the T point, 
the degenerate VB is composed of p x and p y orbitals, 
while the CB is formed by the hybridization of s and p z 
orbitals. 

In Fig(ljc) we also present the electronic band dis- 
persion taking into account spin-orbit coupling (SOC). 
Though clearly the inclusion of spin-orbit interaction 
does not result in a visible change in band dispersion, 




FIG. 3: (Color online) Side and top view for characteristic 
adsorption geometries for alkali, alkaline earth and transition 
metal atoms. 

a band gap of A — 1.5 meV appears at the K point. Due 
to its buckled structure, silicene has a larger SO-induced 
gap than graphene which is of the order of 10 -3 meVP^ 
The calculated band structure and band gap opening (A) 
with S OC is in good agreement with recently reported 

studiespmm 

In Table [T] we also compare structural, electronic and 
vibrational properties of graphene and silicene. It ap- 
pears that in contrast to the general trend that the larger 
the atomic radius, the smaller the workfunction, silicene's 
workfunction is 4.77 eV, while for graphene it is 4.49 eV. 
The calculated values of in plane stiffness^ and cohesive 
energy indicates that silicene is a stable but less stiffer 
material as compared to graphene. Similar to graphene 's 
out-of-plane optical (ZO) phonon mode at 900 cm' 1 , sil- 
icene has a ZO mode at 150 cm -1 . Eigenfrequency of 
the LO and TO modes are degenerate at the T symme- 
try point are found to be 580 cmT x which is almost three 
times smaller than graphene's LO (and TO) modes. 



IV. RESULTS: ADSORPTION OF METAL 
ATOMS 

As shown in Fig. [2] regarding the interaction of silicene 
surface with adsorbates four different adsorption posi- 
tions can be considered i.e., above the center of hexag- 
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FIG. 4: (Color online) Diffusion pathways and barriers of (a) 
alkali atoms and (b) 3d, 4d and 5d transition metals through 
top, bridge, valley and hollow sites. 



onal silicon rings (Hollow site), on top of the upper sil- 
icon atoms (Top-site) , on top of the lower silicon atoms 
(Valley-site), on top of the Si-Si bond (Bridge site). Con- 
sidering the monolayer hexagonal lattice structure of sil- 
icene, it is reasonable to expect the relaxation of foreign 
atoms to one of these adsorption sites. 



A. Bonding Geometry and Migration Barriers 

We first investigate the adsorption characteristics of 
alkali adatoms Li, Na and K on silicene. The alkali met- 
als are highly reactive metals and their chemical activity 
increases from Li to Fr. Characteristic bonding geom- 
etry of alkali atoms are depicted in Figj3] Upon full 
geometry optimization, all alkali atoms Li, Na and K 
atoms favors bonding on the hollow site of the silicene 
layer. Adsorption of alkali atoms does not yield any sig- 
nificant distortion or stress on the silicene lattice. The 
valley site on the low-lying silicon atoms is the next fa- 
vorable site. Though the top and valley site adsorptions 
are also possible, bridge site adsorption of alkalis is not 



possible on a silicene lattice. Therefore the bridge site 
adsorption is a kind of transition state between top and 
valley sites. Structural and electronic properties of alkali 
metals adsorbed on silicene layer are also presented in 
Table [IT] Here the height of adatoms is calculated as the 
difference between the average coordinates of neighbor- 
ing Si atoms and the adsorbate. The distance between 
the adatom and silicene surface monotonically increases 
with increasing atomic size. However, fully conforming 
to graphenepB there is no clear trend in the adsorption 
energies. While adsorption energies of Li, Na and K on 
graphene were calculated to be 1.1, 0.5 and 0.8 eV, re- 
spectively, their binding to silicene lattice is more than 
twice stronger, i.e. 2.4, 1.9 and 2.1 eV, respectively. The 
nature of the alkali-silicene bond will be discussed in the 
next subsection. 

Possible diffusion pathways of the adsorbate atoms on 
silicene lattice can also be deduced from the energy bar- 
rier between hollow, valley, top and bridge sites. In Fig. 
[4] we present the energetics of different adsorption sites. 
For alkali atoms, shown in Fig. |4ja) it appears that the 
most likely migration path between subsequent hollow 
sites passes thorough the nearest valley sites. It is also 
seen that the energy difference between hollow and val- 
ley sites becomes smaller for larger atoms and therefore 
the diffusion of larger alkali atoms are relatively easier. 
Though alkali atoms strongly bind to the silicene sur- 
face, at high temperatures they may diffuse along hollow 
and valley sites when they overcome the energy barrier 
of 140-280 meV. 

Alkaline earth metals are the elements of the periodic 
table having two valence electrons in their outermost or- 
bital. Compared to alkalis, alkaline earth metals have 
smaller atomic size, higher melting point, higher ioniza- 
tion energy and larger effective charge. Strong interac- 
tion between alkaline earths and silicon surfaces is well- 
known and have been used for various silicon etching and 
surface engineering techniques. Therefore one can expect 
strong bonding of alkaline earth's to a monolayer silicon 
surface. Resulting atomic geometries for Be, Mg and Ca 
adsorbed on a silicene sheet are depicted in Fig. [3j Unlike 
alkalis, the hollow site is not the most favorable adsorp- 
tion site for alkaline earth metals. Among these, while Be 
and Mg favor adsorption to a valley site, the Ca adsor- 
bate that has a quite large atomic size (empirically ~ 1.8 
A) prefers bridge site adsorption. Adsorption energy of 
alkaline earth's is slightly higher for alkali metals (except 
Mg). Similar to the Na adsorbate a second row element, 
adsorption energy has a sharp decrease for the second 
row alkali element Mg. As shown in Fig. Qb), the Be 
atom has to overcome a quite large energy barrier (AE 
> 1 eV) for migration from valley to other adsorption 
sites. Sudden increases in migration barrier stem from 
the stretching of silicene lattice by adsorbate which is 
freed only in one direction perpendicular to the surface. 
At high temperatures, migration of Ca atoms through 
bridge and hollow sites may take place by overcoming 
the energy barrier of 125 meV. 
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TABLE II: Calculated values for adatom adsorption on 
(6x6x1) silicene; relaxation sites hollow (H), bridge (B), val- 
ley (V) or top (T), adatom height (h), adsorption energy of 
adatom (-EaiJs); total magnetic moment of the system (fitot) 
in units of Bohr magneton (pb), energy band gap (E g ), dipole 
moment (p), Bader charge transfered from adatom to silicene 
(Pad) an d the workfunction of the optimized structure (<&). 
Metallic and half-metallic structures are denoted as m and 
hm, respectively. 
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We next investigate the adsorption characteristics of 
eight elements of the 3d, Ad and 5d transition metal 
adatoms: Ti, V, Cr, Mn, Fe, Co, Mo and W. Though 
the outermost s orbitals of the transition metals are 
completely filled, because of their partially filled inner 
d orbitals diverse adsorption characteristics for different 
atoms can be expected. Due to relatively small atomic 
radius of all transition metals and having more electrons 
that can participate in the chemical bonding, we can ex- 
pect stronger binding to the silicene lattice. From Fig. 
[3] it appears that while the adsorption of alkalis do not 
cause any significant change in the silicene lattice, transi- 
tion metals are more likely to disturb the nearest silicon 
atoms. 

Bonding of Ti adatom with 4.89 eV to silicene occurs 
with a significant lattice distortion whereas the adsorp- 
tion of Ti occurs on the hollow site of graphene with- 
out disturbing the planar lattice structure. As a conse- 
quence of bridge-site adsorption of Ti, it binds six near- 
est Si atoms strongly by pushing the underlying two Si 
atoms downwards. In order not to exclude possible va- 
cancy formation and adatom-induced fracturing in such 
Ti-adsorbed silicene lattice we also examine the stabil- 
ity of whole structure through molecular dynamics (MD) 
calculations. Ab initio MD calculations show that the 
adsorbed Ti atom remains bounded and neighboring Si- 
Si bonds are not broken after ~1 ps at 300 K. Similar 
to Ti, V adatom is adsorbed on bridge site with 4.32 eV 
adsorption energy. 



Among the 3d transition metal adatoms considered 
here only Cr, Mn, Fe and Co similar to alkalis, are ad- 
sorbed on the hollow site. However, differing from alkalis, 
transition metals have quite strong binding (3.20, 3.48, 
4.79 and 5.61 for Cr, Mn, Fe and Co, respectively) with 
three uppermost Si atoms. Therefore, instead of being 
adsorbed on top a hollow site like alkalis, transition met- 
als are almost confined in the silicene plane. Energetics of 
transition metal atoms on the most favorable adsorption 
sites are also listed in Table [ill 

Note that not only the electronic occupancy but also 
the atomic radius of the atom is important in determining 
the final geometry of the adatom on silicene. Since for 
d orbitals having more than half-occupancy the atomic 
radii are significantly decreased by increasing the number 
of electrons, starting from Cr all 3d transition metals 
prefer to be relaxed on hollow site. To see how the atomic 
radius of transition metal affects the adsorption geometry 
we also perform calculations for other Group 6 elements 
Mo and W. According to the most recent measurement 
of Cordero et alW^ covalent atomic radii of Cr, Mn and 
W are 1.39, 1.54 and 1.62 A, respectively. In this column 
of the Periodic Table, while Cr is adsorbed on hollow 
site, bridge site adsorption becomes more preferable for 
Mo and W due to their larger atomic radii. Effect of 
the atomic radii can also be seen even in the same row 
elements: adatoms Ti and V, that have covalent radius 
larger than that of Cr, are relaxed to the bridge site. 
Thus we can infer that only the transition metal adatoms 
having covalent atomic radii larger than ^1.50 A favor 
bridge site adsorption. 

It appears from Fig. [4^c) that the energy barrier be- 
tween different adsorption sites is relatively large for tran- 
sition metal atoms. The diffusion barrier between the 
most favorable site and the less favorable one is ^0.6 
eV. The most likely diffusion path for Ti, V, Mo and W 
passes from bridge to hollow sites. However, diffusion of 
Mn, Cr, Fe and Co atoms from one hollow site to another 
one may occur via valley sites. 

Furthermore, the effect of spin-orbit interaction on 
the optimized adsorbate-silicene geometry is examined 
for adsorption on (6x6x1) silicene supercell. Compared 
with the optimized geometries obtained excluding SOC, 
SO-induced change in adatom-silicon bond length is just 
on the order of 0.001 A. It is also found that the migra- 
tion path profile and the most favorable adsorption site 
do not change when SOC is included. Therefore, due 
to negligible effect of spin-orbit interaction on structural 
properties, in the rest of our study LDA calculations will 
be employed. 



B. Electronic Structure 

In this section we present spin-polarized electronic 
band dispersion and density of states (DOS) for adatoms 
adsorbed on the most favorable site on the silicene sur- 
face. Since the alkaline earth and transition metal ad- 
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FIG. 5: (Color online) Electronic band structures for perfect, 
Li, Na and K adsorbed silicene. Fermi level is set to zero. 
Occupied bands are filled with yellow color. 




FIG. 6: (Color online) Partial DOS for Be, Mg and Ca ad- 
sorbates and nearest silicon atoms. Fermi level is set to zero. 
DOS is broadened by Gaussian smearing with 0.2 eV. 



sorption significantly disturb the hexagonal lattice sym- 
metry, electronic band dispersion along the high symme- 
try points (r-M-K-r) of perfect silicene may not repre- 
sent the real electronic properties of the whole structure. 
Therefore, a density of states plot covering large number 
of k-points in the BZ is more convenient for a reliable 
description of the electronic structure. 

In Fig. [5j we present electronic band dispersions of al- 
kali metal adsorbed silicene. For the sake of comparison, 
the band structure of the bare silicene is also included. 
Here it is worth to note that the Dirac cone that is lo- 
cated at 2/3 of first BZ is shifted from K to T high sym- 
metry point due to band folding of the (6x6x1) silicene 
supercell. Similar shift can be obtained by band folding 
operations for all (3nx3nxl) supercells (n=l, 2,...). As 
a result of the adsorption of an alkali atom, semimetallic 
silicene becomes metallic due to the donation of ~0.8 e 
charge from the alkali atom into the silicene conduction 
band. Such an attraction of adatom through the hollow 
site of the silicene surface with a significant charge trans- 
fer resembles ionic bonding. Moreover, since the band 
dispersion of perfect silicene is negligibly disturbed by 
the alkali adatoms, it seems reasonable to assume that 
the bonding is of ionic type. Here, in accordance to 
the charge transfer, the Fermi level (Ep) is shifted by 
~0.2 eV. It is seen that all alkali atom adsorption bands 
formed by the hybridization of adsorbate-s states with 
the silicene-p states appear in the vicinity of 0.4 eV. Ad- 
ditionally, a small gap opening in the p bands below the 
Fermi level, which is 6, 13 and 31 meV for Li, Na and K, 
respectively, occur at the crossing point. 

While each isolated alkali atom has a net initial mag- 
netic moment of 1 /is, up-down spin degeneracy is not 
broken upon 0.8 e charge transfer and therefore all the 
silicene+alkali structures are nonmagnetic metals. How- 
ever, as a result of the large charge transfer between the 
alkali atom and silicene, remarkable dipole moment per- 



pendicular to the silicene surface is induced. Calculated 
value of the dipole moment directed from silicene to the 
adatom is 0.30, 0.60 and 0.94 e-A for Li, Na and K, re- 
spectively. Though the amount of the charge transfered 
is almost the same for all alkalis, the induced dipole mo- 
ment is different due to the difference in distance be- 
tween the adsorbate and the silicene layer. As a con- 
sequence of adatom-induced surface charge density re- 
distribution, the workfunction of the silicene surface is 
decreased. Therefore, the trend of the increase in the 
dipole moment follows the decrease in workfunction. 

Alkaline earth metals are highly reactive elements that 
tend to form various compounds by loosing two out- 
ermost shell electrons. Natural silicate and carbonate 
forms of alkaline earths are well-known. As shown in Ta- 
ble [TTJ Be, Mg and Ca adsorption yields a gap opening of 
0.39, 0.48 and 0.17 eV, respectively. It is also seen from 
the partial density of states shown in Fig. [6] that the top 
of the valence band (VB) of Be adsorbed silicene is due 
to the hybridization of s and p z orbitals of the nearest 
silicon atoms and the alkaline earth metals. However, 
at top of the VB main contribution comes from Mg-s, 
Si-s and Si-p z hybridization, whereas p z orbital of Mg 
does not mix with surrounding silicon states. The bot- 
tom of the conduction band (CB) of both Be and Mg is 
formed by hybridization of p xy states of adsorbates with 
p xy and p z states of silicene. Differing from Be and Mg, 
adsorption of Ca on silicene does not yield appearance of 
Ca states around ~Ep- While the states at the VB edge 
arise mainly from the silicene p z orbitals, the CB edge is 
composed of s, p xy and p z orbitals of Si. 

Since atomic d states of transition metal atoms have 
comparable energy values to that of their s valence states, 
the d shell of these atoms is partially filled. The net mag- 
netic moment of isolated transition metals are nonzero 
and behave like small magnets unless the d shell is com- 
pletely filled. Though the d shell electrons are located 
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FIG. 7: (Color online) Partial DOS for Ti, V, Cr and Mn 
adsorbates and nearest silicon atoms. Fermi level is set to 
zero. DOS is broadened by Gaussian smearing with 0.2 eV. 



close to the nucleus like the characteristic core electrons, 
they can spread out much further like valence electrons. 
Therefore, transition metals with their partially filled d 
shells are relaxed to different sites on silicene and we can 
expect diverse electronic properties upon their adsorp- 
tion. We depict the spin-polarized electronic density of 
states of those transition metal adsorbed silicene in Fig 

m 

Ti adsorption on silicene has a quite different charac- 
teristics as compared to other transition metals. Upon 
the adsorption of a Ti atom that has an initial magnetic 
moment of 4.0 fis in its isolated state, on bridge site of 
silicene the structure becomes a ferromagnetic metal with 
2.0 [Ib net magnetic moment per supercell. Since we use 
a quite large supercell such that the adatom-adatom in- 
teractions are negligible, describing the adatom levels in 
terms of orbitals, instead of bands, is more appropriate. 
As a result of the strong bonding, the silicene lattice is 
significantly distorted, Ti atom binds to six nearest Si 
atoms and the degeneracies of the Ti-c? states are com- 
pletely broken. Metallic bands of the structure originates 



only from the Ti-d 2 .2_ y 2 (f) orbitals and silicene becomes 
half metal when Ti adatom is adsorbed. 

For the case of V adsorption on silicene, the struc- 
ture becomes semiconductor with 2.7 \ib net magnetic 
moment per supercell. It is seen from Fig [7] that while 
only Y-d XVtyz xz {\) states have a significant contribution 
to VB maximum, bottom of the CB is derived from Si- 
p{\.) and ~V-d X y^y Z ^ xz (\r). The contribution of the other 
eigenstates at the vicinity of F,p are enhanced just be- 
cause of the smearing of the DOS plot. Cr adsorption 
also turns the semimetallic silicene into a half metallic 
material with 4 fj,B per supercell. The up spin (minor- 
ity) bands of the Cr adsorbed silicene are semiconducting 
with a direct band gap of ~0.4 eV, whereas the down spin 
(majority) band shows metallic behavior. The bands in 
the vicinity of are mainly derived from the hybridiza- 
tion of Gv-d X y tVZyXZ {\.) and Cr-d x 2_ y 2 Q.) orbitals and the 
Si-p(4_) states. However, differing from Cr, the adsorp- 
tion of Mo and W that occurs on bridge site results in a 
nonmagnetic ground state. While Mo-adsorbed silicene 
becomes metallic due to Mo-d xy , yZyXZ states at E^, W- 
adsorbed silicene is a nonmagnetic semiconductor with 
bandgap of 0.02 eV. Note that the half-metallic behavior 
for Ti- and Cr-decorated silicene can be quite important 
for potential use in spintronics. 

The adsorption of Mn atom in close proximity to the 
hollow site of buckled hexagonal lattice results in semi- 
conducting behavior with a bandgap of 0.24 eV. The de- 
generacy of f and J. spin states are broken due to the exis- 
tence of a net magnetic moment of 3.0 \ib- Since Mn, Fe 
and Co atoms have the smallest atomic size of the consid- 
ered metal atoms here, they are the most closely bonded 
one among all hollow site adsorbed atoms. Both the ad- 
sorption of Fe and Co result in metallic ground state, 
that originates from the d xy ^ VZtXZ {\,) and d z 2(\.) states at 
E^, with 2 and 1 fig net magnetic moment, respectively. 
As a consequence of the small adsorption height, the in- 
duced dipole moments arc much smaller than those of 
the alkali adatom adsorbed silicene. 

The adsorbate-induced magnetic moments induced in 
the silicene layer lead to an exchange-splitting especially 
in d eigenstates. The calculated value of the exchange- 
splitting is 0.4, 2.0, 2.1, 1.9, 0.5 and 0.4 eV for Ti, V, 
Cr, Mn, Fe and Co, respectively. However, bridge-site 
adsorbed Mo and W adatoms exhibits no splitting due 
to their nonmagnetic ground state. 

To reveal the correlations between adsorption energy 
(F<Ads), workfunction ($) and induced dipole moment (p) 
we also present "J-E^s and p-$ plots as shown in Fig [8] 
Note that the adsorption of large transition metals that 
occurs with large binding energy changes the silicene 's 
workfunction negligibly. Especially for 3d transition met- 
als having more than half filled rf-subshell, the larger the 
adsorption energy the higher the workfunction. However, 
it appears that even with a small coverage (1/72) of sil- 
icene surface by alkali metals, a significant decrease of 
the work function can be produced. For alkalis and alka- 
line earths, there is almost a linear correlation between 
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FIG. 8: (Color online) (a) Plot of the adsorption energy E^ds 
versus the workfunction <!> and (b) the workfunction <!> versus 
dipole moment p. 



F,Ads and $. Because of the ionic nature of the alkali- 
silicene bonding, the workfunction linearly depends on 
the atomic size and therefore one can expect a signifi- 
cant decrease (>1 eV) in the workfunction for adsorption 
of larger alkalis on silicene. It is also seen that, while 
multivalent adatoms do not follow a particular trend in 
dipole moment, for alkali atoms workfunction shift ex- 
hibits a linear decrease with increasing dipole moment. 



CONCLUSIONS 



the diversity of adsorption geometries and its high surface 
reactivity. As a consequence of its buckled lattice struc- 
ture, all the metal atoms bind strongly to silicene. While 
the binding energy of alkali and alkaline earth metals 
is 1-3 eV, transition metals have larger binding energies 
of 3-7 eV. Our diffusion path analysis shows that com- 
pared to graphene migration of metal atoms on a silicene 
lattice is more difficult and requires to overcome higher 
energy barriers. On the other hand, diffusion of metal 
atoms Li, Na, K, Mg and Ca may occur at moderate tem- 
peratures. Depending on the adatom type, semimetallic 
silicene can show metal, half-metal and semiconducting 
behavior. It was also noted that as a result of the high 
surface reactivity the direction of the charge transfer is 
always from metal adsorbate to silicene. However, the 
existence of charge donation and the resulting adatom- 
induced dipole modify the workfunction of silicene con- 
siderably. Especially the linear decrease in $ for larger 
alkali adatoms is promising for tunable enhancement of 
electron and ion emission which is appealing to silicene- 
based thermionic converters and cathodes. Our findings 
also suggest that the half-metallic ferromagnetic nature 
of Ti- and Cr-decorated silicene have a great potential 
for silicon-based spintronic device applications. Moreover 
the existence of a tunable bandgap opening in silicene by 
alkaline earth metal adatoms is highly desirable for its 
use in nanoscale optoelectronic device applications. 

The extension of our investigation to a detailed study 
to include the effect of spin-orbit coupling and of the 
substrate interactions on the adsorption characteristics 
of silicene is planned for future studies. 



In this study, we presented a first-principles investiga- 
tion of the adsorption characteristics of alkali (Li, Na and 
K), alkaline earth (Be, Mg and Ca) and transition metals 
(Ti, V, Cr, Mn, Fe, Co, Mo and W) on silicene. Inves- 
tigation of the interaction of silicene with metal atoms 
has a significant importance because of its fundamental 
relevance to applications in catalysis, batteries and na- 
noelectronics. 

We found that silicene has quite diff erent adsorption 
characteristics as compared to graphen d 65 * 66 ! because of 
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